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Abstract 

Two kinds of asymmetry arise from the interference of the Born amplitude and the box-type am- 
plitude corresponding to two virtual photons exchange, namely charge-odd and one spin asymme- 
tries. In case of unpolarized particles the charge-odd correlation is calculated. It can be measured 
in combination of electron muon and positron muon scattering experiments. The forward-backward 
asymmetry is the corresponding quantity which can be measured for the crossed processes. In the 
case of polarized muon the one-spin asymmetry for annihilation and scattering channels has been 
calculated. The additional structure function arising from the interference is shown to suffer from 
infrared divergencies. The background due to electroweak interaction is discussed. 

PACS numbers: 



1 



I. INTRODUCTION 



The motivation of this work is to give an accurate description of the process ep, — > 
epi^y) ee — > pp{"f) in frame of QED, in order to provide a basis for the comparison with 
experimental data. High precision experiments on the processes ee — > rf and ee — > pp are 
planned in future c — r facilities. Moreover, the possibility of colliding e/x beam facilities has 
been discussed in framewok of programs on verification of Standard Model (SM) prediction. 

The obtained results can also be applied, as a realistic model, to electron (positron) 
scattering on a point-like hadron (proton), based on arguments given in this paper. 

It is known that in Born approximation the differential cross-section of elastic proton- 
electron scattering 

e(pi) + p(p) -> e(pi) + p(p') (1) 

can be expressed in terms of two proton form factors, F 12 (g 2 ), which are functions of a single 
argument, the momentum transfer squared, q 2 = t. 

Taking into account two (and more) photon exchanges, (TPE), leads to a generalization 
of the Born picture, namely the amplitude of ep scattering depends on two Mandelstam 
variables, the total energy s and t. The virtual photon Compton scattering amplitude is a 
rather complex object, which can be expressed in terms of 12 chiral amplitudes. Nevertheless, 
taking into account parity conservation and omitting the terms of order m e /m M , (which are 
responsible for chirality violation) reduce the number of relevant amplitudes to three l|: 

M w = ^L^hXpi) x Hp') [fiM^ - 1 ^^f-iA + ~^(M)(pi + #)(p+ A*]"(p), 

with q = pi — p[, s = 2pxp. 

The explicit calculation, given in this paper, permits one to extract the individual con- 
tributions Fx, F2, and F3, in frame of QED. The infrared divergency is cancelled when the 
relevant soft photon emission is correctly taken into account. 

Charge-odd and backward-forward asymmetries appear naturally from the interference 
of one and two photon exchange amplitudes in frame of QED and SM due to Zo-boson 
exchange in Born approximation. But at the energy range reachable at c — r factories, the 
relevant contribution of SM type is |2j : 

z - r a v a a w 5 ■ 1(T 5 , 3 < < 5 GeV. (2) 
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which is quite small compared to QED effects. 

The accuracy of results given below is determined by 
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and the contribution of higher orders of QED a /it « 0.5%. Moreover we assume that all the 
velocities of the final heavy particles are finite in the annihilation as well as in the scattering 
channels. This is the reason why Coulomb factors are neglected. 

Our paper is organized as follows. In sections |H] the annihilation channel and the scat- 
tering channel ee — > are considered in sections HP and ITTT1 respectively. In section HVl we 
take into account the soft photon emission and construct charge-odd and forward-backward 
asymmetries. In section we analyze the crossing relation between two channels. Explicit 
form for additional structure G% for annihilation channel is given in section IVTl In section IVm 
the one-spin asymmetries are investigated. The results are summarized in the Conclusions. 

II. PROCESS e+ + e- ->p+ + M -( 7 ) 

At first, we consider the process of creation of fi + n~ pairs in electron-positron annihila- 
tion: 

e + (p+)+e-(p^) ^ fi + (q+)+fi-(q-). (4) 
The cross-section in the Born approximation, can be written as: 

with s = (p_|_ + p-) 2 = AE 2 } (3 2 = 1 — ^j-, E IS the electron beam energy in center of mass 

reference frame (implied for this process below), m, m e are the masses of muon and electron, 

c = cos#, and 9 is the angle of /i_-meson emission to the electron beam direction. 

The interference of the Born amplitude 

iAttq: 
s 

with the box-type amplitude Mg, results in parity violating contributions to the differential 
cross section, i. e. the ones, changing the sign at 9 — ► ix — 9,. As a consequence of charge-odd 
correlations we can construct: 




FIG. 1: Feynman diagrams for two-photon exchange in ee — * fifi process: box diagram (a) and 
crossed box diagram (b). 

Here we take into account as well the emission of an additional soft real photon with energy 
not exceeding some small value AE, so that A(9, AE) is free from the infrared singularities. 

Part of the results presented here were previously derived by one of us (E. A. K.) in Ref. 
0|, and partially published in 4|. 

There are two box- type Feynman amplitudes (Fig. [TJ. We calculate only one of them, 
the uncrossed diagram (Fig. with matrix element 



M a = ia 2 J 



d ku(q-)Tv(q + ) x v(p+)Zu(p-) 



(A)(Q)(P + )(P_) 

(A) = (k - A) 2 - ml (Q) = ( k -Q) 2 -m 2 , (P±) = (k T P? - A 2 , (7) 

with A-" photon" mass and 

T = ^ a {k -Q + m)nfp, Z = jp(k - A)j a , 

A = -(p + -p_), Q = ±(q + -q_), P = ±(p ++ p_). (8) 

We will assume 

m 2 = -(1 - f3 2 ) ~ s t~-u. (9) 

The explicit form of kinematical variables used below is: 

A 2 = -P 2 = - S -, Q 2 = -I s /3 2 , a= AQ =-(«-*), 
u = (p_ - q + ? = -|(1 + P 2 + 2/3c), t = (p_ - g_) 2 = -£(1 + /3 2 - 2/3c). (10) 

The contribution to the cross section of the amplitude arising from the crossed Feynman 
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diagram (Fig. Hh), M&, can be obtained from M a by the crossing relation 

da a (s,t) da b (s,u) 



dfi M df^ 



(11) 



which has the form 

da a (s,t) Pa 3 



40, 2^ R ^- < 12 > 



with 



f d 4 k 11 1 

R ( s > *) = J ^ (A)(Q)(p + )(p_) 4 Tr + m ) T ^+ ~ m ) ?aJ x 4 Tr (P+ Z P-^) ■ (13) 

The scalar, vector and tensor loop momentum integrals are defined as: 

j. j . j _ f 1' h^ky , . 

' ^ ~ 7 m 2 (A)(Q)(P+)(P_) 1 J 

Using symmetry properties, the vector and tensor integrals can be written as: 

J„= Ja- A m + J q -Q^ (15) 

= Kog^ + K P P^P U + K Q Q»Q V + K A A^A U + K X {Q^A V + Q V A^. (16) 
The quantity R(s, t) can be expressed as a function of polynomials Pj as: 

R = P\J + P 2 Ja + P 3 Jq + P 4 K + P 5 K A + P 6 K Q + P 7 K P + P S K X , (17) 

where the explicit form of polynomials is given in Appendix A. Using the explicit expression 
for the coefficients Ja, K x (See Appendix B) we obtain 

R(s, t) = 4(cr - A 2 )(2a - m 2 )F + lQ(a - A 2 )(a 2 + (A 2 ) 2 - m 2 A 2 )J 

+4[(A 2 ) 2 - 3AV + 2a 2 - m 2 a]F Q + 4[2(A 2 ) 2 - 2A 2 ct + 2a 2 - m 2 A 2 ]F A 
+4[(A 2 ) 2 + A 2 a + m 2 A 2 ]G Q + 4[-(A 2 ) 2 + a 2 - 2m 2 A 2 }H Q , (18) 

with the quantities F 4- Hq given in Appendix B. Finally the charge-odd part of differential 
cross section has the form 

fda% rt (s,t)\ _ a 3 (3_ vann 



dfi M ) odd 27rS 
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Li 2 (z) = - / — ln(l-x) 



(20) 



is the Spence function. The quantity D ann — T>y in suffers from infrared divergences, which 
will be compensated taking into account the soft photons contribution (see below). 

III. SCATTERING CHANNEL 



Let us consider now the elastic electron muon scattering 

e(pi) + M -> e(Pi) + Kp') 

which is the crossed process of (j3J). The Born cross section is the same for the scattering 
of electrons and positrons on the same target. Taking the experimental data from the 
scattering of electron and positron on the same target (muon or proton), one can measure 
the difference of the corresponding cross-sections which is sensitive to the interference of the 
one and two photon exchange amplitudes. For the case of proton target, in the Laboratory 
(Lab) frame, the differential cross section as a function of the energy of the initial electron, 
E and of the electron scattering angle, 9 e , was derived in Ref. Q]: 



da ep 



a 2 cos 2 ^ 1 
AE 2 sin 4 \ p 



F 2 4- tF 2 



E ^ M +2rF 2 ,tan 2 " e 



1 + r 



2E 



p = 1 H sin — 



-t 



M 

E 2 



2 



• 2 
2 , S111 ~ 



(21) 



m 2 ' Am 2 m 2 p 2 

and it is known as the Rosenbluth formula. The Sachs electric and magnetic proton form 
factors, F E and F M are related to the Pauli and Dirac form factors by F E — F x — tF 2 , 



Fm = F\ + F 2 . For the scattering on muon, one replaces F% = 1, F2 = and Eq. 
becomes 



da% a 2 (s 2 + u 2 + 2tm 2 ) 



s = 2pip = 2mE, t = —2pip[, u = —2pp[ 



dVl 2m 2 p 2 t 2 p 

The charge-odd contribution to the cross section of e/i— elastic scattering is: 



dVt e 
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In 2 - + 2Li 2 1 + — 
i V m 
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with the help of the following relation: 

1 



+ 2m 2 - 



su 



In 



ln^ 



m? + s m? + u 



m z F c 



4,/r(l+r) 



7T 



+ ln(4r) lnx + Li 2 (-2 



tx) — Li 5 
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where 



x 



y/1 + T + y/7 
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IV. SOFT PHOTON EMISSION 



In this section the emission of soft real photons in the Lab reference frame for eu- scat- 

n 

tering is calculated. Following Ref. the odd part of cross section 



da soft 



a 



• 2 



d 3 k ( p\_ _ pi_\ ( j/_ p_\ 



(26) 



da An 2 J u \p[k p t kj \p'k pk J 5oj w<Ae 

must be calculated in the special reference frame So, where the sum of the three- momenta 
of proton and of the recoil proton is zero z = k + p = 0. Really, in this frame, the on-mass 
shell condition of the scattered muon S[(z — k) 2 — m 2 }, z = p\ + p — p[ does not depend on 
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FIG. 2: H(s,cos#) for E = 5m (dashed line) and E = 10m, m is muon mass. 



the direction of the emitted photon. The photon energy can be determined as the difference 
of the energy of the scattered electron and the corresponding value for the elastic case: the 
maximum value of the photon energy Ae in the Sq frame is related with the energy of the 
scattered electron, detected in the Lab frame AE as (see [6], Appendix C), 



Ae = pAE. 



(27) 



The calculation of the soft photon integral with u < Ae can be performed using t'Hooft and 
M. Veltman approach (see Section 7). We find 

a 3 (s 2 + u 2 + 2tm 2 ) 



da soft 



2 In p In 



soft 



!) + 2Li 2 (l-^),„=A 
px J \ xj —u 



7r 2m 2 p 2 t 2 
-2Li 2 (l - 

which is in agreement with Ref. 

The sum (da™ r * + dcr|°^) odd has the form 

. a 3 (s 2 + u 2 + 2tm 2 

odd 



(2pAEf 
X 2 x 



+ V 



soft ( ) 



'da virt 
dtt P 



dag^ 

dn P 



2irm 2 p 2 
Re 



t 2 



-2 In p In 



(2pAEf 
—tx 



f 2 T) sc 



T) sc 

s 2 + u 2 + 2tm 2 soft 



and it is independent from the photon mass A. 

The function 5 is shown in Fig. |2]as a function of cos^g for given E/m . 
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The ratio between the difference and the sum (corresponding to the Born cross section) 
of the cross sections for e^fi scattering is: 



foe H~*e 0(7) _ foe+ fj,-^e+ n(j) 
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( J (T e+/i-ve+/i(7) _j_ foe- n-*e- /rift) 

The odd contributions to the differential cross section for the process e + + e" 
due to soft photon emission, has the form: 

d 3 A; 



(30) 



(da 



e + e — >/^ + /^ (7)\ 
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I odd 
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a: 



4vr 2 



P- P+ 



q+ g_ 



p_A; p + A; I \q+k q-k 



• (31) 
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Again, the integration must be performed in the special frame S°, where p+ + P- — q+ 
q_ + k = 0. In this frame we have 



(q- + k) 2 - m 2 = 2(£_ + u)u « 2mu = (p + + p. 
E-e, 



m 2 = AE(E - 



m . 
— As. 
2E 



In the elastic case E — e+ = and the photon energy in the Lab system is 

m 



AE = e 



el 



e. 



2E 



:Ae. 



The t'Hooft-Veltman procedure for soft photon emission contribution leads to: 
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(34) 



(l + /3)(l-/3c), 
/ 2/3(1 + c) 



;i+/3)(i+/3 C )J- (35) 

The total contribution (virtual and soft) is free from infrared singularities and has the form 
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The quantity $(s,cos0) is presented in Fig. El 

The relevant asymmetry can be constructed from © 
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FIG. 3: <E>(s,cos0), for s = Wm 2 (dashed line) and s = 20m 2 , m is muon mass. 



V. CROSSING SYMMETRY 



In this section we formally consider the relations between the kinematical variables in 
the scattering and in the annihilation channel, e + + e~ — > p + p. The reduced form of the 
differential elastic ep scattering cross section, commonly used, is defined as a re( i = rFlj+eF^ 
and it is related to the differential cross section by: 

do- a 2 cos 2 % 1 1 



— — — aM&red, CM ~ 

ail 



4E 2 sin 4 |pe(l + r 



1 + 2(1 +t) tan 2 f ' 



(38) 



where e is the polarization of the virtual photon, and varies from e — 0, for 9 e — ir to e — 1, 
for e e = 0. 

The crossing relation between the scattering channel e + p — > e + p and the annihilation 
channel e + + e _ — > p + p consists in replacing the variables of the scattering channel s = 
2pip = 2Em and Q 2 = —t according to 

m 2 + s -> t = -2£ 2 (1 - /5c), Q 2 -> -s = -4E 2 , cos0 = cosp_V-, (39) 

where is the angle of the antiproton with respect to the incident electron, in the center of 
mass system (CMS). 

The following relation holds, for the annihilation channel: 

(t - it) 2 



cos* 6 



s(s-AM 2 ) 

On the other hand, in the scattering channel, one has: 



, s + t + u = 2m . 



1 + e cot 



2 e e 



l-e 



1 + T 



+ 1 



[s - u 



Q 2 (Q 2 + AM 2 ) 
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, Q 2 



s + u. 



(40) 



(41) 



Therefore one proves the validity of the crossing relation: 



ooBfl = yi±|=y, (42) 

based on the analytical continuation from the annihilation channel to the scattering one. 
This relation was derived in Refs. (s[ 0|. Using this relation and the property of 
the 27 contribution to the annihilation cross-section ^j^(^)^ = — (j^C^ — i- e - 

( = cos6'/(cos 2 6, s) the authors of Ref. [I| built the ansatz for 27 contribution to 
ep-elastic scattering 

^(e-p -> e-p) = yf(y 2 , Q 2 ); f(y 2 , Q 2 ) = c (Q 2 ) + y 2 Cl (Q 2 ) + y%(Q 2 ) + .... (43) 

This property follows from the change of the sign of the contribution for virtual and real 
photon emission when the (s <-» u) transformation is applied (see Eqs. ()24|28I29|) and 
relation (@TJ)). 

This form of the contribution of the I7 ® 27 interference to the differential cross section 
derives explicitly from C-invariance and crossing symmetry of electromagnetic interactions 
and excludes any linear function of e for a possible parameterization of such contribution. 

Let us note that not only the elastic channel must be taken into account: the interference 
of the amplitudes corresponding to the emission of a photon by electron and by proton must 
be considered, too. 

Evidently, the relations derived above are valid for the considered processes with electrons 
and heavy lepton setting respectively F E = F M = 1. 



VI. DERIVATION OF THE ADDITIONAL STRUCTURE: ANNIHILATION 
CHANNEL 

Let us start from the following form of the matrix element for the process e + (p + ) + 
e (p ) ~ * A i+ ( ( ?+) + ^ {l-) i n presence of 27 exchange: 

M 2 = — v(p + )^u(p-) x u(q-) (G xlll - + A-G.AQ,) v(q + ), (44) 

S \ 7Yt S J 

where the amplitudes Gi are complex functions of the two kinematical variables s, and t. 

To calculate the structure G 3 from the 27 amplitude (see Eq. (J7|)), both Feynman dia- 
grams (Figs. and^i) must be taken into account. Similarly to Section II, only one of 
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them can be calculated explicitly (the uncrossed one), whereas the other can be obtained 
from this one by appropriate replacements. 

To extract the structure G 3 we multiply Eq. (|44|) subsequently by 



u{p-)l\v{p+) x v(q + )j x u(q- 
u(p-)Qv(p+) x v(q+)u(qJ), 
u{p-)Qv(p + ) x v(q + )hu(q_ 

and perform the summation on fermions spin states. 
Solving the algebraical set of equations we find 
1 



/3 4 sin 4 fl 
° 2 = /? 4 sin 4 # 
° 3 = p 3 sin 4 d 



1 - (3 2 ) {A a (3 sin 2 9 - 8C a cos 9) + AB C 



/3 2 (l + cos 2 #)]} 



A a (3 2 sin 2 9 cos 9 - 8B a cos 9 + A(3C a (l + cos 2 9) 



with 



B a 



d 4 k 
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d 4 fc 



1 



Tfl A ] 

^ (A)(Q)(P + )(P_) ^ rr( ^' g) X 4 Tr + m)T( ^ 



d 4 A; 



1 



^ (A)(g)(p + )(p_) ^ rr( ^-^ } x \ Tr ^ + m)T( ^+ - m)A 



The explicit value for is: 

3 /? 3 (l-c 2 ) 2 



{^ Q (l-c 2 )/3 3 (l-/3c) 



+ ^H Q (3 2 (l - c 2 ) [c(-3 + 5(3 2 ) - /? - (3c 2 



2 ) [c(-3 + 5/3 2 ; 

2 i o,q4 , „2 q2Io a a2\ r>a„n n 



(45) 



{ (8P a + A a (3 2 sin 2 0) (1 - (3 2 cos 2 0) - AC a (3 cos [2 - (3 2 ( 1 + cos 2 6*)] } 



(46) 



m)] , (47) 



(48) 



+ F A c [1 - A(3 2 + 2(3 A + c 2 (3 2 (3 - A(3 Z ) - 2(3c(l - 2(3' 
+ F Q (3 [-c 2 + (3c {-\ - A(3 2 c 2 + \c 2 ) + (3 2 {-\ 
- 2Js/? 2 c(l-c 2 )(l-/3 2 )(l-/?c) 

+ Fc [l + (3 2 c 2 - 2/? 4 - 4/3 4 c 2 + (3c(-3 + A(3 2 + 2(3* + (3 2 c 2 )] } . 
The contributions from the crossed Feynman diagram can be obtained from Eqs. EfSlbv: 

(A\B\C b ) crossed = -[A a ,B a ,C a (cos9 - cos 9)] uncrossed . (49) 
As one can see, in the quantities G\, G 2 , and G 3 infrared divergencies are present. 
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VII. ONE SPIN ASYMMETRY 



Let us consider now the process of electron interaction with a heavy lepton (point-like 
proton). For clearness the expressions are written for the proton case. For the interaction 
of electrons with heavy leptons, \i or r, one should take Ge = Gm = 1 and use the relevant 
mass replacement. 

The target spin asymmetry for heavy fermions production process e + (p + ) + e~(p_) — ► 
p(q+) +p(q-) (in CMS frame) is defined as 

dcr T — da^ 

d^W = ( 5 °) 

where a is the proton polarization vector, n = (<f_ x p_) / \q_ x p_\ is the unit vector nor- 
mal to the scattering plane, d<r^ is the cross section of processes with proton polarization 
vector a, da^ is the cross section of processes with proton polarization vector —a. Thus the 
denominator of the left hand side in Eq. ()50|) is the unpolarized cross section of process 
e + e~ — > p + p which is well-known jlO| (for the case of proton-antiproton creation): 



dVL 4s 



(l + cos 2 6) \G M \ 2 + (l - P 2 ) \G E \ 2 sin 2 e] , (51) 

with (3 = \Jl — ^p- is the velocity of proton in cm. frame, s is the total energy square and 
9 is the angle between vectors q_ and p_. 

The difference of cross sections in (I50j) is originated by the s-channel discontinuity of 
interference of the Born-amplitude with TPE amplitude 

da T - da 1 ~ Re Utiastic ■ Atpe + A dastic ■ A+ PE ) . (52) 



Using the density matrix of final proton u(p)u(p) = (p + M)(l — 75a) one gets 



, (Anay (2m) 



■\2 



Re T. (Static ' a tpe + A aaMc ■ Af PE ) = 32i i-A -Re{Y), 



Y 



x-Tr 
4 



(A)(Q)(+)(- 
(p - M) (- T5 a) 7a (p' + M) 7, (fc - Q + M) lll 



Performing the loop-momenta integration the right hand side of Eq. (|53|) can be expressed 
in terms of basic integrals (see Appendix B) 

Re(Y) = 4M(o, A, Q, P)Im (F Q - G Q + H Q ) , (54) 
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FIG. 4: Asymmetry R n for the case of structureless proton for energies s = 5 GeV 2 (dashed line) 
and s = 15 GeV 2 (solid line). 

where (a, A,Q,P) = e^ upo a^A u Q p P a = ( v /i/2) 3 (an)/?sin 6. Using the expressions listed in 
Appendix B we have: 

lm { F Q - G Q + H Q ) = tm = JL (i^ In i±| - 2) . (55) 

Thus, after standard algebra, the following expression for spin asymmetry can be obtained 
for the processes e + + e~ — > p + p: 

M /#(/?) sin 

^ = 2 "7i 2-/3 2 sin^ (56) 
and it is shown in Fig. 01 as a function of 9 at several values of s for the case of structureless 
proton. 

Such considerations apply to the scattering channel when the initial protons is polarized. 
Similarly to (fHKJ) one finds 

Im.(F Q -G Q + H Q ) = -j^. (57) 

(note that the s-channel imaginary part vanishes for the crossed photon diagram amplitude). 
The contribution of the polarization vector appears in the same combination 

1 ME 2 
(a,A,Q,P) = -{a,p,p h q) = —^-sm6(an). (58) 
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FIG. 5: Asymmetry T n in case of structureless proton for E = 5 GeV (dashed line) and E = 10 
GeV (solid line). 

The single spin asymmetry for the process e~ + p — > e~ + p (the initial proton is polarized) 
has the form: 

da 1 * — da^ 



do"T -)- dcr^ 



(an)T n (59) 



with 



r » = 2^JT^( 1 + T »^ " s ta 4 (60) 

This quantity is shown in Fig. Elfor the case of structureless proton as a function of 9, for two 
values of s. The asymmetry decreases when the c.m.s. energy growth, so on experiment it is 
useful to measure the asymmetry near the treshold of proton (or heavy lepton) production. 



VIII. CONCLUSIONS 

We calculated QED radiative corrections to the differential cross-section of the processes 
e + + — > fi + + ^(7), e ± + fi — > e ± + ^(7), arising from the interference between the Born 
and the box-type amplitudes. The relevant part of soft photon emission contribution, which 
eliminates the infrared singularities, was also considered. 
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According to the previous discussion, the box Feynman diagram, calculated in QED 
frame, can also be applied to lepton-hadron interaction. In the present work, we consid- 
ered structureless Born approximation for annihilation and scattering channels. A realistic 
case with Born and box diagram depending on form-factors will be considered further (in 
preparation). 

Angular asymmetry, charge asymmetry as well as target spin asymmetry were calcu- 
lated. These quantities are free from infrared and electron mass singularities. Numerical 
applications show that these observables are large enough to be measured (see Figs. I2JEJ). 

The charge-asymmetry properties of radiative corrections in the annihilation channel 
induce non-trivial terms in the cross section due to crossing symmetry relations (see (|43j) ) . 

The parametrization (J44j) for the contribution to the matrix element arising from box- 
type diagrams in terms of three additional functions Gi(s, t), i = 1, 2, 3 suffers from infrared 
divergencies, which can be eliminated by taking into account soft photon emission expressed 
in terms of structures G±, G2, G3. This procedure results in replacing ln(m/A) with In A. 

The results obtained here, for the processes e^y — » e ± p(7) are par ticularly interesting in 
view of the experiments planned at Novosibirsk |lll| and at JLab^|l2| as well as e + + e" 
NN("f), which can be investigated at Frascati jl3j and Bejing 

IX. APPENDIX A: TRACE CALCULATION. 

The explicit expressions for the polynomials Pj are: 



0. 



Pi = 


8{-(A 2 ) 3 


- AV + 2a 3 + [(A 2 ) 2 - 2A 2 a]m 2 }, 


P2 = 


16[(A 2 )V 


- a 3 + AVm 2 ], 


P3 = 


8({2(A 2 ) 2 


a - 2a 3 + m 2 [(A 2 ) 2 + 2AV - a 2 } + A 2 m 4 }, 


Pa = 


8[5(A 2 ) 2 - 


-6A 2 (T + 5(T 2 -5A 2 m 2 ], 


P 5 = 


8[(A 2 ) 3 - 


2(A 2 )V + AV-m 2 A 2 )], 


Pe = 


8{(A 2 ) 3 - 


2(A 2 )V + AV + m 2 [-(A 2 ) 2 - 2AV + 2a 2 ] - 2A 2 m' 


Pi = 


8[-(A 2 ) 3 - 


■f 2(A 2 )V- AV + (A 2 ) 2 m 2 ], 


P 8 = 


8{-(A 2 ) 3 


+ (A 2 )V - 3AV + 3a 3 - m 2 [(A 2 ) 2 + 3A 2 a]}. 



(61) 
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X. APPENDIX B: USEFUL INTEGRALS. 



In the calculation of ep scattering we use the following set of scalar integrals with three 
and four denominators 



—i 



F A -- — 
Fq 



ir 2 J (A)(P+)(P_ 
d 4 k 



1 



— + - ln z — 

2 mi 



7T 

1 

» = -% 



(Q)(P+)(P- 

2 



2 2 



d 4 fc 



(A)(g)(p 4 

■ 2 -t 



+ Li 2 



/ , m — t , m 
+ 2 In — + ln — 

1 



In 



7T Z 



m 
A 2 



'1 + /3 N 



d 4 fc 



- Li 2 



(a)(q)(p_; 



In 2 



-sJ-G 
2 



2(m 2 - t) 



2 ki- 



rn — t 



2Li 2 



1 

2(m 2 - 1) 
t 



In 2 



+ In — In — 



mi 



m c 



In 2 



m 



m 



t 1 , 9 m 2 
2 - + -ln 2 — + 2Li 2 

2 mf 



J = — 



d 4 £; 



2 In 



m z — t 



m 



+ In— ] In — 



tt 2 7 (A)(Q)(P + )(P_) s(™ 2 -*) 
The terms proportional to m 2 /s, m 2 /m 2 were neglected. Notations follow Ref. 
The vector integrals with three denominators are: 



1 

in 
H A 



k^drk 



(a)(q)(p + 
i / 



H P P» + H A A» + H Q Q", H Q = -\n 



m 2 — t 



m c 



m 2 — t 
H + - 



2 

mi 



m 2 + t , m 2 — t' 



t 



In 



m 2 — t 
k^d 4 k 



mJ (A)(P + )(P_) 
fcM 4 jfe 



' m 2 , m 2 — A 
ln — + 2 ln — , 

= GaA m , G A = if-21n4 + iln 2 4 + ^ 

s \ mi I mi d , 



G 



1 



mi (Q)(P+)(P_) ' ~ Q s-4m 2 

Four denominator vector and tensor integrals were defined in (I14J) . The relevant coefficients 
are: 



-2 ln — + sF<; 
m z 



(62) 



(63) 



Ja " 2d 
1 

Jq = 2rf 



(P + P Q )a- A 2 (P + P A ) 



F = \sJ-G, d = A 2 Q 2 



a 
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1 

2~a 



K A 

K P 

K Q = 
K x = 

where we used 



a{F - G + Hp + H A + H Q ) + H A {a - A 2 ) - H Q {a - Q 2 ) 
+2P 2 J A (A 2 - 2a) + A 2 G A - Q 2 G Q - 2p 2 Q 2 J Q 

-7T1 \Q 2 o{G -F-H p -3H a + 6P 2 J a ) 
2aa L 

+(A 2 g 2 + a 2 )(H A - 2P 2 J A - G A ) - Q\H Q - 2P 2 J Q - G c 



1 



2a(H A - 2P 2 J A + H p + l -F- l -G) + Q 2 (H Q - 2P 2 J Q - G Q ) 



2P 2 a 

-A 2 (H A -2P 2 J A - G 



2ad 



- A z aA P + 2A 4 A A + (o- 2 - 2A 2 Q 2 )A Q 



^(aA P + Q 2 A Q -2A 2 A A ), 



(64) 



A A = H A + 2A 2 J A -G A , Aq = Hq + 2A 2 Jq — Gq, A p = F -G+ H P + 3H A + 6A 2 J A . (65) 
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